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ABSTRACT 

Vertical magnetic fields have been known to exist in the internetwork region 
for decades, while the properties of horizontal magnetic fields have recently been 
extensively investigated with Hinode. Vertical and horizontal magnetic fields 
in the internetwork region are considered to be separate entities and have thus 
far not been investigated in a unified way. Wc discover clear positional asso- 
ciation between the vertical and horizontal magnetic fields in the internetwork 
region with Hinode. Essentially all of the horizontal magnetic patches are as- 
sociated with the vertical magnetic patches. Alternatively, half of the vertical 
magnetic patches accommodate the horizontal magnetic patches. These horizon- 
tal patches are located around the borders of the vertical patches. The intrinsic 
magnetic field strength as obtained with the Stokes V line ratio inside the hor- 
izontal patches is weak, and is in sub-equipartition field regime {B < 700 G), 
while the field strength outside the horizontal patches ranges from weak to strong 
(kG) fields. Vertical magnetic patches are known to be concentrated on meso- 
granular and supergranular boundaries, while the horizontal magnetic patches 
are found only on the mesogranular boundaries. These observations provide us 
with new information on the origin of the vertical and horizontal internetwork 
magnetic fields, in a unified way. We conjecture that internetwork magnetic fields 
are provided by emergence of small-scale fiux tubes with bipolar footpoints, and 
the vertical magnetic fields of the footpoints are intensified to kG fields due to 
convective collapse. Resultant strong vertical fields are advected by the super- 
granular fiow, and eventually form the network fields. 

Subject headings: magnetic fields — convection — Sun: photosphere — Sun: 
granulation 
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Introduction 



In the quiet regions of the Sun, strong vertical fields are concentrated around the 



boundaries of the supergranulation, and are referred to as networ^ 



i fields. The intrinsic 



magnetic field strength is determined to be 1-2 kG (IStenfio 



19731 ). and is comparable 



to the field strength in the plage regions. The re g ion in side the network is referred to 



as the internetwork region. 



Livingston fc Harvey 



( 1l97ll ) found that the internetwork 



region has weak magnetic fields. Since then, spectropolarime tric observations have been 



carried out to explore the internetwork magnetic 



Lin fc Rimmele 



1999 



Socas-Navarro fc Lites 



ields (e.g. 



Lin 



1995 



Lites et al 



1996 



2004 ). Most of the observations have used 



only a ci rcular polarization signa l, and thus, the magnetic field inclination could not be 



deduced. 



Khomenko et al. 



(120031 ) carried out the full spectropolarimetric observation using 



19951 ) than the often used 



Fe I 1565 nm lines, which have higher Zeeman sensitivity (iLinl 
Fe I 630.2 nm lines, and found that the magnetic fields in the internetwork region have a 
broad range of inclinations with peak essentially at vertical. 



le.g, 



Martin 



Schriiver et al.l (Il997( ) concluded that all of the fiux in the network is generated locally 



19901), and only the small net fiux density is the result of dispersal from 
decaying active regions. They also point out that the fiux supphed by the emergence of the 
ephemeral regions is the most important source for the quiet-Sun network and that they 
cannot put constraints on the amount of fiux that is injected on scales smaller than that of 
the ephemeral regions. Even with these pioneering works, the properties and the origin of 
the internetwork fields have been far from understood, and high-resolution observations of 
the internetwork region is critically important to clarify the generation mechanism of the 
quiet Sun magnetic fields. 



The Solar Optica 


Telescope spectropolarimeter (SO''. 


VSP; 


Tsuneta et al. 


2008b 


Suematsu et al. 


2008 




Ichimoto et al. 


2008 




Shimizu et al. 


2008 


) on board Hinode 
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( iKosugi et al.ll2007( ) has allowed us to measure the full Stokes spectra with high polarization 
sensitivity and with high spatial resolution. The SOT/SP revealed the ubiquitous 
horizontally inclined magnetic fields with the intrinsic field s trength of a few hundred 



G in the internetwork regions (jOrozco Suarez et al. 



2007a 



Lites et al. 



internetwork hori zontal magnetic fields were originally reported by 



Lites et al. 



2008). Such 



fll996h and 



De Pontieul ( 12002| ). Owing to the high spatial resolution and the high polarimetric sensitivity 



of SOT/SP, tl 


le pro 


jerties of horizonta 


fields have been studied in greater detail (e.g., 


Centeno et al. 


2007 




[shikawa et al. 


2008 


Jin et al. 


2009; 


Martinez Gonzalez & Bellot Rubio 


2009 




Ishikawa et al. 


2010 


). The horizontal magnetic fields in the internetwork region 



were also confirmed and inve s tigated ( 



mission ( iBarthol et al. 



2010 



Danilovic et al 



Solanki et al. 



2010; 



20101) with the recent S unrise /IMaX 



Martinez Pillet et al. 



2010|). The 



horizontal magnetic fields a re highly transient with the lifetime raging from 1 to 10 minutes 



( llshikawa fc Tsuneta 



2009al ). The se horizontal magneti c fields are observed in the quiet Sun 



as we. 



2008ah . T 



1 as in a weak plage region ( llshikawa et al 



20081) and a polar region (ITsuneta et al. 



lis ubiquity is consi stent with that of seething horizontal magnetic fields reported 



earlier by 



Harvey et al 



( 120071 ) using the SOLIS and GONG instruments. The discovery of 
the ubiquitous horizontal fields provides us with new insight to better understand the quiet 
Sun magnetism in more general perspective. 



An important question on the properties of the internetwork 



concerned with their inclination distributions. 



Orozco Suarez et al 



magnet ic fie 



(l2007bl ) and 



ds is 



Lites et al. 



(120081 ) report ed that the ho r izonta l ma gnetic fields are domina nt in the internetwork 



region, while 



Beck fc Rezaei 



(120091) and 



Khomenko et al. 



( 120031 ) found the dominance of 



the vertical magnetic fields. Some authors pointed out that magnetic field vectors are 



rando mly oriented and there is no pre 



region ( iMartmez Gonzalez et al. 



2008 



erence of a specific inclination in the internetwork 



Asensio Ramos 



2009 



Bommier et al 



20091 ). Most of 



these authors exploit the probability density function (PDF) of the inclination angles and 
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the field strengths. Such PDFs represent one aspect of the internetwork magnetic fields. 
We, however, stress that the original spectropolarimetric images have rich information 
which cannot be represented by simple PDF. For instance, one can address, based on these 
images, whether the horizontal field patches are somehow associated with the vertical fields 
or are a completely different entity from the vertical ones. In this paper, we study the 
spatial distributions of the vertical and horizontal magnetic fields in an internetwork region, 
and aim to address a question on how these horizontal and vertical fields are spatially 
distributed and possibly related. 

We emphasize the importance of the spatial distribution of the vertical and the 
horizontal magnetic fields in terms of the convective patterns with different scale sizes 



to understand the re 



ationship of the magnetic fields with the convective motion. 



Ishikawa fc Tsunetal ( 120101 ) found that the horizontal fields appear with clear cellular 
pattern of a typical scale of 5-10" and that the cellular structure coincides in position 
with the negative divergence of the horizontal flow field, i.e., mesogranular boundaries 
with downfiows. As to the internetwork vertical fields, similar concentrations have been 



observed ( iDommguez Cerdena 



20031 ). Mesogranular flows would play an important role on 
the generation and distributions of the internetwork magnetic fields. A natural question 
is whether there are similar concentrations of the internetwork fields on the supergranular 
scale. In this paper, we present the spatial distributions of the vertical and the horizontal 
magnetic fields in terms of mesogranular and supergranular convection. 

One of the major issues about the inter network fields is the ir origin. One possibility 
is the "recycling" mechanism pointed out by 



Ploner et al. 



(120011 ). The remnants of active 



region flux are continually recycled by granular and supergranular convection to form 



the network as well as internetwork fie 
process as suggested by 



ds. Another possibility is that local dynamo 



Cattaned (119991 ) works for the production and maintenance of the 
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i nternetwork fields. An "exp losion" of a rising flux tube from deeper convection zones 



( iMoreno-Insertis et al. 



19951) may generate weak magnetic fields in the upper convection 



zone that would be the source of the quiet-Sun magnetic fields. 

The field strength of the horizontal fields in the internetwork region is a few 
hundred G, which is below the equipartition field strength corresponding to the granular 
convective motion. The histograms (PDF) of the intrinsic magnetic field strengths of the 
horizontal fields are exactly the same in the polar region, the quiet sun, and a weak plage 
region in spite of the signif i cant difference in the amount of the vertical magnetic flux 



(llshikawa &: Tsuneta 



2009a : 



Ito et al 



20101 ). This ubiquity combined with the unique PDF 
and their sub-equipartition field strength implies that the horizontal magnetic fields are 
generated by the local dynamo process, in which the granular convective motion amplifies 
the magnetic fields by folding and stretching. By simultaneously focusing on both horizontal 
and vertical fields, we try to reveal the origin of the internetwork magnetic fields in a unified 



way in this paper. 



2. Observations and Data Reduction 

2.1. Hinode Observations 

A quiet region with size of 302" x 162" at the disk center was observed using the 
Solar Optical Telescope Spectropolarimeter (SOT/SP) on March 10, 2007. The SOT/SP 
provides us with the full Stokes /, Q, U and V spectra covering the Fe I 630.15 nm (Lande 
factor (76302 = 1.66 7) and Fe I 63 0.25 nm ((^eaoi = 2.50) lines with the wavelength sampling 



of 2.14 pm (llshikawa et al. 



2010l ). Each slit position has an exposure time of 4.8 s. The 



noise level (Tq estimated from the continuum in the Stokes Q and U spectra is 1.2 x 10 '^J, 



and that in the Stokes V spectrum is 1.1 x 10 Jc for a single wavelength pixel, where I, 
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is the continuum intensity. The scanning step is O'.'IS and the spatial samphng along the 
slit is &!16. The dark current, the flat field corrections, and the polarimetric calibration 



were performed using "sp_prep" software available in SolarSo: 



extensively usee 


by 


Asensio Ramos 


2009 



many authors (e.g., 



Viticchie et al 



Orozco Suarez et al. 



't. The identical data set was 



2007a 



Lites et al. 



2008 



20101 ) to explore the quiet Sun magnetism. In this 
paper we analyze a 67'.'5 x 72" subfield, which covers one supergranular cell, and contains 
both network and internetwork regions (Figure [1]) . Approximate scan time of the subfield 



is 36 min, which is much shorter than the lifetime of the supergranular cell. 



2.2. Selection of valid pixels 

The polarization signals in the internetwork region are typically small as compared 
with those in the network regions and active regions, and they are likely to be contaminated 
by the photon noise. It is important to carefully choose pixels for analysis that are 
minimally affected by the noise. First of all, we calculate the wavelength-integrated degree 
of polarization (Qtot, ^tot and Vtot) defined as, 

_ /^^; Q{X)dX 

/J d\ 



Utot ^ -^Hxr^, (2) 



_ I,y{X)dX A>(AMA 
i.j^yx - I.J-dX ' 

where and Xr indicate the wavelength limits of the integration in blue and red wings 
of the Fe I 630.25 nm line, respectively, and Aq represents the wavelength of the line 
center. The wavelength Aq corresponds to the minimum value of Stokes / for Qtot and 
Utot, and to the zero-crossing wavelength of the Stokes V profile for V^ot- Here we take 
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Fig. 1. — A 67! 5 x 72" sub-field selected for analysis, (a) continuum intensity map (b) map 
of net linear polarization Ltot as defined by equation (jlj). (c) map of circular polarization 
Vtot as defined by equation ([3]). 
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Ab = Ao — 17.1 pm, and Xb = Xq + 17.1 pm. In figure [21 the histograms of Qtot and Vtot in 
our analyzed subfield (figure [T]) are shown with dashed lines. Since the histograms of Utot 
are the same as those of Qtot, we do not show them. 

To estimate the noise for Qtot, fAot, and Vtot we obtain the wavelength-integration of 
the Stokes Q, U, and V spectra {Q[^^, Ul^^, and V^'^J in the wavelength range minimally 
affected by line polarization (i.e., continuum range). The integrations are carried out over 
the same wavelength width of 36.4 pm. Then, we derive the noise level o"tot which is the 
standard deviation of Q[ot, t/tot? ^ind V^^^, and the value is estimated to be 2.6 x 10~^Ic for 
Qtot, 2.5 X 10~^ Jc for f/tot, and 2.7 x 10~^/c for V^ot- Such wavelength-integrations reduce 
the noise level, and the sensitivity to weak magnetic field significantly increases. The 
histograms of Q^t) Kot corresponding to the noise distributions are shown with solid 
lines in the top panels of figure O Six vertical lines in each upper panel indicate 1, 2, and 
ScTtot- For our analysis, we choose the pixels with Qtot or Utot or Vtot above corresponding 
3atot to avoid the pixels with false polarization signal (i.e., noise). 



Viticchie et al 



(I2OIOI ) report that in the internetwork region, there are pixels which 
have heavily asymmetric Stokes profiles such as one lobe and three-lobe Stokes V profiles 
in the Hinode data. On e example is the footpoint regions of the horizontal magnetic fields 



m 



Ishikawa et al. 



( I2OIOI ). Even if the maximum amplitudes of the Stokes profiles are large 



compared with the noise level uo defined and estimated in Section [2?T| wavelength- integrated 
polarization signals are small, and may not surpass the threshold of 3crtot- We include 
these pixels with high amplitude and low-integrated polarization signal in our analysis. 
We measure the maximum amplitudes of each Stokes profile of the Fe I 630.25 nm line 
(Qmax, Umax and Vmax) to find the pixels with sufficient polarization signal induced by the 
magnetic fields. In the bottom two panels of figure El the histograms of Qmax and \4iax in 
our analyzed region are shown with dashed lines, while the solid lines show the histograms 
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of maximum amplitudes of Stokes Q and V spectra in the wavelength range without any 
polarization signal (i.e., continuum range). These solid lines indicate the noise distributions 
of Qmax and Vmax- Eight vertical lines represent l-4cro for each Stokes profile. Pixels with 
Qniax and Vmax abovc 4(To are essentially free from noise. Since the histograms of ?7max are 
the same as those of Qmax, we do not show them. 

In summary, we choose the valid pixels with Qtot, ^tot, or Vtot larger than Sutot or with 
Qmax, f^max, or \4iax larger than Aa^. To further reduce the effect of noise, we choose patches 
that contain at least 4 contiguous pixels, all of which meet the above criteria. The size of 4 
pixels corresponds approximately to the instrument's spatial resolution. According to this 
selection procedure, 10.1% of the subfield has sufficient linear polarization signals (Stokes 
Q and/or t/), and we refer to these pixels as horizontal field pixels (hereafter referred to 
as HF pixels). Likewise, 43.8% has sufficient circular polarization signal (Stokes V), and 
these pixels are referred to as vertical field pixels (VF pixels). Some pixels are identified as 
both the VF and the HF pixels. In total, 46.9% of the subfields are valid pixels in terms of 



above se^ 
35.5% of 



ection criteria. This number (46.9%) for Stokes Q, U o t V is much larger than 



Orozco Suarez et al. 



(12007br ) and 29% of IViticchie et all fl2010f ). They investigated 
the same data set, but analyzed pixels only with the maximum amplitudes larger than 4.5 
times their noise levels. 

In the subfields, we identify 1085 HF patches and 1793 VF patches, respectively. 
Among 1793 VF patches, 901 VF patches have positive polarity, and 892 VF patches have 
negative polarity. In figure E] (a), the VF patches with negative and positive polarities are 
shown in black and white, respectively. The red contours in the figure indicate the HF 
patches. 
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Fig. 2. — Upper panels: Histograms (broken lines) of wavelength-integrated polarization 
signals (Qtot and Vtot) and histograms (solid lines) of the noise (Qtot Kot)- Lower panels: 
Histograms (broken lines) of the peak polarization signals (Qmax and Knax), and histograms 
(sohd lines) of the peak noise levels. 
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2.3. Conversion of polarization signal to magnetic flux density 



As a result of measurement noise, inversion of Stokes profiles for internetwork fields may 
not converge, thus compromising the uniqueness and accuracy of the solution. Furthermore, 
the asymmetric profiles often seen in the internetwork region make the proper inversion 
more difficult. Indeed, Milne-Eddington inversion, which has been widely used for the 
analysis of the SOT/SP data, cannot deal with the heavily asymmetric Stokes profiles. 
More po^ yerful inversion techniques such as SIR (Stokes Inversion based on Response 



function, 



Ruiz Cobo &: del Toro Iniesta 



Atmosphere, 



Sanchez Almeida 



19921 ) and MISMA (Micro-Structured Magnetic 



19971 ) are useful to interpret the Stokes profiles with large 
asymmetries. However, in these techniques the number of fitting parameters increases, and 
we may potentially face the severe problem of the uniqueness of the solution. 

Therefore, in this paper, we do not use the inversion technique to retrieve the 
properties of the internetwork magnetic fields. Instead, we follow a procedure to convert 
the wavelength-integration of the polarization signal in to the longitudin al and transverse 
magnetic flux density. Following previous studies (e.g., iLites et al.ll2008l ). we refer to these 
two measurements as "apparent flux density", B^^^ and -B^p, and use the units of Mx cm~^ 
for them, while the unit of gauss is used for intrinsic fleld strength. This procedure does not 
allow us to resolve flne magnetic structure within the resolution element, which is usually 
represented by the magnetic flUing factor /. and B^ are average fleld strengths over 



the resolution element, and are also interpreted as the magnetic fleld strength assuming 
that / = 1. 

The wavelength-integrated circular polarization Vtot deflned by equation is converted 



to B^ . We also deflne the net linear polarization Ltot which will be converted to B^ 



app> 



(4) 
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Here At and are defined in Section [2^2] 



Lites et al.l (Il999[ ) and iJin et al.l ( 120091 ) employed the inversion results for pixels with 
stronger polarization signal primarily from network regions in the quiet Sun to calibrate 
the relationship of wavelength-integrated polarization signals to the apparent flux densities. 
However, the intrinsic field strengths of the network fields are typically stronger than 
~1 kG, and the calibration curve derived from the network fields may not be applied to the 
internetwork regions which harbor hG magnetic fields (e.g., 



Orozco Suarez et al. 



2007a 



bD 



The approach to determine the 
and is the following (e.g.. 



longitudinal and transverse magnetic flux density 



Lites et al. 



20081 ). We calculate the synthetic Stokes 



profiles based on the Milne-Eddington atmosphere. For the calculation, the intrinsic field 
strength B is set to be smaller than 1 kG, the stray light fraction, a (a = 1 — /) to be 
smaller than 0.6, and the absolute value of Doppler velocity wlos to be less than 2 km s~^. 
Here we choose the field strengths, with which the weak field approximation is valid. We 
employ relatively large filling factor (/ > 0.4) since the small filling factor with weak field 
strength does not produce enough Stokes signals. Many synthetic profiles were generated 
for < < 1 kG, 0.4 </ < 1.0, -2.0 < wlos < 2.0 km s"^ and < 7 < 180°, where 7 is 
the inclination angle with respect to the line-of-sight (LOS) direction. Following equations 
([3]) and (jlj), we calculate Vtot and Ltot from the synthetic profiles. The longitudinal and 



transverse flux density B^^^ and i?^^ are obtained by 



-app 



B 



app 



fB cos 7, 



and 



^Jpp = y^Bsin-f. 



(5) 



(6) 



The wavelength-integrated degree of polarization Vtot and Ltot are compared with B^^^ and 



-^app to derive the calibration curves optimized for the internetwork fields (Figure |3]). Since 



the observing region is located at disk center, and i?fpp correspond to the horizontal 
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and vertical magnetic flux density, respectively. 

The measurement noise (Tq at each wavelength position is not canceled out, and 
coherently contributes to the integration over the wavelength for Ltot (equation @). 
Figure |3] (c) shows the noise distribution of the net linear polarization signal Ltot, which is 
derived with equation but with the integration carried out over the continuum range. 
We set 5j,p = for the mean value Ltot = 1-51 x 10 "^/c, which is indicated by the solid line 
in figure |3] (c). In general, the distribution function p{l) with I = ^/ q"^ + u^, where q and u 
have the Gaussian distributions with the same standard deviation a {p{q) = exp( — ^) 
and p{u) = exp( — ^)) is written as p{l) = ^exp(— ^). Thus, the mean value of / is 
slightly deviated from the peak of the distribution as shown in figure [3] (c). Note that this 
calibration will underestimate the magnetic flux density where the intrinsic field strengths 
surpass kG strengths, as it does in network regions. 

Figure H] shows the spatial distributions of longitudinal (vertical) and transverse 
(horizontal) magnetic flux density, -Bfpp and -B^p having been obtained through the 
procedure described above. The scales for both components are saturated at 250 Mx cm~^ 
in order to compare the two maps in the internetwork region. 

We identify the pixels with Vtot > 0.5% and the area >90 pix^ as network pixels, 
which are marked with yellow contours in figure [6] (b). HF patches that overlap with the 
network pixels are indicated with green contours in figure [6] (b). These HF patches may be 
the part of the canopies of the network fields. We do not include these HF patches in the 
subsequent analysis. The red contours in figure M (b) indicate the horizontal magnetic fields 
(HF patches) in the internetwork regions. 
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3. Close association of horizontal magnetic patches with vertical magnetic 

patches 

In this section, we examine the spatial relationship between the horizontal and vertical 
magnetic fields. Figure M (a) and figure [7] (a) indicate that both horizontal and vertical 
magnetic compon ents together form the cellular structures with a spatial scale of 5" — 10" 



(ILites et al. 



20081 ). We also notice that the HF patches overlap with the VF patches. 



As mentioned in section [2l2l the total number of VF pixels in total is 4.3 times as large 
as that of the HF pixels, and 69.4% of the HF pixels overlap the VF pixels. Alternatively, 
48.7% of the total number of the VF patches has cospatial HF patches. The size of the 
VF patches is apparently larger than that of the HF patches (figure E] (a)). This would be 
because the sensitivity of circular polarization to detect VF patches is much higher than 
that of linear polarization to detect HF patches. Thus we detect only the core portion of 
the HF patches. Even though many HF patches overlap with the VF patches, they are not 
contained in the VF patches, and tend to be located around the edge of the VF patches 
(figure [6] (a)). We quantitatively examine where th e HF patches a. r e loca ted with respect to 



the VF patches by means of the procedure used in llshikawa et al. 



fl2007h 



The area density D{L) of the HF patches in the internetwork region at the distance 
L {L < inward, and L > outward) from the borders of the VF patches with width dL 
(figure [5]) is given by 

where A^aii(-^) is the total number of pixels that form the VF patches between L and 
L + dL, and A^hf(-^) is the total number of pixels occupied by the HF patches located in 
the distance between L and L + dL. Here we take dL = 1 pixel (0'.'16). Note that Ngxi[L) 
and A^hf(-^^) are not obtained for a single VF patch, but for all the VF patches. 



- 16 - 



The procedure to obtain Ni{L) {i = all, HF) with the outward direction (L > 0) is as 
follows. Let us consider the case with one VF patch shown in black in figure O (a). A^aii(-^) 
is the total number of pixels in the gray area of figure E] (b) and Ai"aii(L) is the total number 
of pixels occupied by the HF patches in the same gray area. First, we obtain the total area 
of the black and white region (figure [5] (b)) Saii{L), and the total area of the HF patches in 
the black and white region S}ip{L). Second, we obtain Si^ii{L + dL), which is the total area 
of the black, white, and the gray region (figured (b)), and Suf{L + dL), which is the total 
area occupied by HFs in the same black, white, and gray region. Finally, the areas Ni{L) 
with i = all, HF are given by Ni{L) = Si{L + dL) — Si{L). Similarly we obtain Ni{L) with 
the inward direction (L < 0). 

Figure [8] (a) shows the area density D{L) of the HFs derived with this procedure. There 
are more HFs around the boundary of the VF pixels; the area density of the HFs reaches 
its maximum point at the inward distance of L = — 0'.'16 (1 pixel) and rapidly decrease over 
L ~0'.'6. This clearly shows a remarkable association of the smaller HF patches with the 
larger VF patches at their borders. 

We detect 1085 HF patches in the internetwork region, and 53% of them have circular 
polarization signals (Kot) of both polarities inside the HF patches and their vicinity 
(within outward distance of one pixel ~0'.'16). The HF patches that have one polarity 
comprise another 43% of the HF patches. The rest of the HF patches (4%) d o not have 



any circular polarization (Kot)- A similar result is reported by 



Danilovic et al 



torn using 



the Sunrise/IMax instrumen t. Horizo n tal m agnetic fields without acco mpanying vertical 



components are reported by 



Jin et al. 



Jin et al. 



fl2009h for 



(120091 ). Such HFs were found by 
~20% of the detectable events, a percentage that is much higher than the one from our 
own analysis (4%). In our analysis, we chose the pixels for which polarization signals are as 



weak as possible, yet still have minimal infiuence from noise. This careful treatment should 
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have enabled us to find a larger nun iber of H 
internetwork region if the results of I Jin et al 



patch es without vertical components in the 



( 120091 ) were correct. In Figure [7] (a), HF 



patches with Vtot of both polarities (bipolar HF patches) are shown with red contours, while 
HF patches with Vtot of one polarity (unipolar HF patches) are shown with green contours. 
The HF patches without Vtot (isolated HF patches) are indicated with blue contours. 



4. Properties of horizontal magnetic patches 

We investigate whether the properties of the HF patches such as their magnetic flux 
density, their size, and their locations depend on the presence of the vertical magnetic 
components. Figure |9] (a) indicates the normalized histograms of the maximum unsigned 
vertical magnetic flux density inside the HF patch. For bipolar HF patches (red lines), the 
histogram includes the maximum value of the positive magnetic flux density and that of the 
unsigned negative magnetic flux density. In the weak end of vertical magnetic flux density 
(<~100 Mx cm^^), there is no significant difference between unipolar HF events (green 
lines) and bipolar ones (red lines). However, more unipolar events harbor relatively strong 
vertical magnetic flux density (> 100 Mx cm~^). These events would represent the canopy 
structures of the vertical magnetic elements, and one example is seen in the unipolar HF 
patch at (27", 46") in figure [7] (a). 

Figure M (b) shows the histogram of the maximum horizontal magnetic flux density 
in the bipolar (red lines), the unipolar (green), and the isolated HF patches (black). Here 
the isolated HF patches refer to those without detectable circular polarization signal. 
The maximum horizontal magnetic flux density of the bipolar HF patches ranges up to 
350 Mx cm~^ with an average value of ~157 Mx cm~^, while the maximum horizontal 
magnetic flux density of the unipolar HF patches ranges to 200 Mx cm~^ with an average 
value of ~112 Mx cm~^ (see also tabled]). On average, the maximum horizontal flux 
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Fig. 3. — Calibration curves derived from tiie syntlietic Stokes profiles based on Milne- 
Eddington atmosphere for weak fields to convert the wavelength-integrated polarization sig- 
nals (T4ot s-nd Ltot) to magnetic flux density for longitudinal magnetic components (a) and 
transverse components (b). (c) Noise distributions of Ltot- 



Table 1: The properties of HF patches. The values in brackets indicate the total number of 
analyzed HF patches. SD represents standard deviation. 









Bipolar 


(574) 


Unipolar 


(467) 


Isolated (43) 


AU (1084) 
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Fig. 4. — Longitudinal {left) and transverse [right) magnetic flux density, -Bfpp and -B^p, 
for the a 67'.'5 x 72" subfield shown in figure [H Color scales for both panels are the same and 
are saturated at 250 Mx cm~^. We set -Bfpp = in pixels which are not identified as the VF 
pixels, and set B'[ = in pixels which are not identified as the HF pixels (Section l2.2p . 



-20- 




Fig. 5. — (a) Schematic diagram of a VF region, (b) Schematic diagram to obtain Ni{L) 
(i = all, HF) with the outward distance L from the edge of the VF region. L is the outward 
distance from the edge of the VF region. See text for details. 
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density of the bipolar HF patches is larger than that of the unipolar HF patches. The flux 
density of the isolated HF patches is essentially below 150 Mx cm~^ with an average value 
of 91 Mx cm~^, and is much weaker than that of HF patches associated with a vertical 
component. 

The histograms of the area of these patches are shown in figure [3 (c). The area of 
the bipolar HF patches distributes over the wide range, while the area of the unipolar HF 
patches is essentially limited to 1.0 arcsec^ and 74% of the unipolar HF patches are smaller 
than 0.2 arcsec^. The area of all the isolated events is smaller than 0.2 arcsec^. The number 
of HF patches generally increases with decreasing size, and about 20 % of the HF patches 
have a size smaller than 0.1 arcsec^ (4 pix), which is comparable to the SOT/SP spatial 
resolution (0'.'32 x 0'.'32). This indicates that there are a number of HF patches that are not 
resolved with SOT/SP. 

Figure M (d) shows the normalized continuum intensity distributions for bipolar and 
unipolar HF patches, and the VF patches (blue line) in the internetwork region. It appears 
that the continuum intensity distributions for both bipolar and unipolar HF patches have 
peaks slightly brighter than that for the whole subfield (dashed line), and are higher 
than the intensity distribution of the whole subfield at intermediate intensities (~0.95 - 



~1.05). This is consistent with the observations ' 


lat t 


[F patches are statistically located 


in the outer boundary of t 


le granules ( 


Lites et al. 


2008 


Centeno et al. 


2007; 


Ishikawa et al. 


2008; 


Ishikawa & Tsuneta 


2009b: 


Danilovic et al. 


201oi 


. We do not display the continuum 



intensity distribution of the isolated HF patches, which is noisy due to the small number 
of pixels. It appears that (1) the continuum intensity distribution of VF patches has a 
peak darker than that of the whole subfield, and is higher at the lower intensities (~0.85 - 
~ 0.95) and that (2) the continuum intensities of VF patches are darker than those of HF 
patches, and VF patches could be located closer to the intergranular lane. 



- 22 - 



The unipolar HF events have lower horizontal magnetic flux densities, and a smaller 
area than the bipolar events. Both bipolar and unipolar events tend to be located in the 
periphery of granules. As to the isolated HF patches, the horizontal magnetic flux densities 
are low and their areas are small compared with the bipolar and unipolar HF patches. 
Although some unipolar HF patches correspond to the canopy structures of the vertical 
magnetic fields, it is conceivable that the some unipolar and isolated HF patches are 
essentially the same phenomena as the bipolar patches; they have the loop structure with 
bipolar footpoints. Since the sensitivity to horizontal field is low compared with sensitivity 
to the vertical field, only cores of the HF patches attached to one of the footpoints are 
detected as unipolar HF patches. This would result in the smaller size of the unipolar HF 
patches. The candidate for such HF patches is found at (25", 43") in figured (a). When 
the HF patches have a size closer to the resolution element of the telescope, their circular 
polarization signals with opposite polarities corresponding to their footpoints start to be 
canceled out. As a result, the residual Stokes V signal becomes weak, and it is possible that 
only circular polarization signals with one dominant polarity are left, or all the Stokes V 
signals are completely canceled out. 

5. Estimation of horizontal and vertical magnetic fluxes 

5.1. Global magnetic flux balance of VF and HF patches 

In this section, we compare the total magnetic fluxes to clarify which component 
dominates the internetwork region, vertical or horizontal. We derive and compare the total 
vertical and horizontal magnetic flux in the pure internetwork region. In figure [6] (c), the 
cross mark indicates the center position of a supergranular cell, and we define the region 
within a distance of 15" from the center (orange circle) as the pure internetwork region. 
Some vertical patches crossing over the orange line are also included in the analysis. We 
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Fig. 6. — (a) VF patches with negative and positive polarities are shown in black and white. 
The red contours represent HF patches, (b) Circular polarization map (Vtot) saturated at 
Vtot = ±0.01. The red contours are the same as the panel (a). Yellow contours indicate the 
network pixels. Green and red contours indicate the HF patches in the network region, and 
the HF patches in the internetwork region. In panels (c) and (d), circular polarization map 
(Vtot) saturated at Vtot = ±0.005 is shown and the red contours indicate the HF patches, (c) 
The cross mark indicates an approximate center position of a supergranular cell. The region 
enclosed by the orange circle is used for the estimation of magnetic flux in the internetwork 
region. Dark blue contours indicate the VF patches associated with the HF patches, while 
light blue ones indicate the VF patches which do not overlap with the HF patches, (d) The 
blue contours indicate the VF regions, while yellow contours indicate the network magnetic 
fields. 



- 24 - 




Fig. 7. — Enlargement of the center part of the internetwork region. The orange circle 
is the same as that in panel (c) of Figure [61 (a) VF patches with negative and positive 
polarities are shown in black and white. Bipolar HF patches are shown with red contours, 
while unipolar ones are shown with green contours. HF patches without vertical magnetic 
components are shown with blue contours, (b) The circular polarization map (Vtot) saturated 
at Vtot = ±0.005 is shown and the red contours indicate the HF patches. 
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Fig. 8. — (a) Area density D{L) [eq.©] of the HF patches in the internetwork region as a 
function of the distance from the boundary of the VF patches. The distance of L = corre- 
sponds to the edge of the VF patches, and the negative and positive distances mean inward 
and outward directions, respectively, (b) Area density of the VF and HF patches in the 
internetwork region as a function of the distance from the center position of a supergranular 
ceU. Sohd hues indicate the running mean with 4 pixels (~ 0'.'6). 
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Fig. 9. — (a) Histograms of maximum vertical magnetic flux density in the HF patches, (b) 
Histograms of maximum horizontal magnetic flux density in the HF patches, (c) Histograms 
of the HF patch area, (d) Histograms of the normalized continuum intensities in the HF 
region. 
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consider the HF patches that are located inside the orange circle and overlap with the VF 
patches (dark blue contours in figure [H] (c)). The VF patches without HFs are shown with 
the light blue contours in figure in](c). 

The total unsigned vertical magnetic fiux within the orange circle is given by 
^ \B cos 'yf\dsdl = ^ \B!^p^\dsdl, where B is the intrinsic magnetic field strength, / the 
filling factor, ds the scanning step size of O'.'IS, and dl is the pixel scale along the slit of 
C/ie. The total fiux is estimated to be 3.0x10^^ Mx cm and the average fiux density 
is 8.3 Mx cm~^. If we consider only the VF patches that have the HF patches, total 
vertical magnetic fiux is estimated to be 2.7x10^^ Mx cm~^, and the average fiux density 
is 7.5 Mx cm~^. We reach the remarkable conclusion that most (~91%) of the vertical 
magnetic fiux are contributed by the VF patches accompanied by the HF patches. The 



mean fl ux density obtained here is comparab le to the one obtained by 



fl2007bl . 9.5 Mx cm'^) and 



Lites et al 



(12008 



Orozco Suarez et al 



11.7 Mx cm-^l 



Following 



Lites et al. 



(120081 ). we also derive the spatial average horizontal flux density 



defined by ^ B^ dsdl/ ^ dS, where ^ dS is the total area enclosed by the orange circle. 



The average horizontal fiux density is 7.1 Mx cm ^ which is comparab 



derived above and much smaller than the one obtained by iLites et al. 



e to the vertical one 
(120081 . 60 Mx cm-2). 



Unlike the VF case, the horizontal magnetic fiux cannot be estimated by -B^p times 



the pixel area, since the magnetic fields are directed perpendicul ar to the LOS 



f all the 



2010h . the 



identified horizontal magnetic patches have a fiux tube structure (llshikawa et al.l 
cross section wdz of the individual fiux tubes instead of the apparent size should be used 
for the estimation of the total magnetic fiux, where w is the lateral size of the fiux tube 
and dz the thickness along the LOS. Assuming that the fiux tube is spatially resolved and 
the magnetic filling factor / is the fraction occupied by the magnetized plasma along the 



LOS in the line-forming layer, the total magnetic fiux of an individual event is estimated 
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by B sin'-fwdz ~ B sinjiufL, where Bsinj is the average horizontal intrinsic field strength 
in the HF patch and L is the length of the line forming region (~300 km). Using equation 



(I6]), the total magnetic fiux is written as B^^^^wL, where i?^p is the average transverse 
(horizontal) magnetic fiux density in the HF patch. Since / is not derived in our analysis, 



we use the average filling f actor f ^ 0.15 of the HFs in the quiet Sun obtained with 



Milne-Eddington inversion fllshikawa fc Tsunete 



2009al ). We fit a spatial ellipse to all the 



HF patches and use the minor axes for the lateral size of the fiux tube w. The total 
magnetic fiux of the HF patches within the orange circle taking into account their cross 
sections is estimated to be 4.7x10^^ Mx. This value is smaller than the total unsigned 
vertical magnetic fiux (3.0x10^^ Mx) by a factor of 6. Since the positive and negative 
magnetic fiux are grossly balanced inside the orange circle, total magnetic fiux for the 
VF patches will be three times larger than that for the HF patches. This value could be 
larger or smaller, since we assume the nominal filling factor of 0.15 and the thickness of the 
line-forming layer of 300 km, and roughly estimate the cross section of the fiux tubes. 



5.2. Local magnetic flux balance of VF and HF patches 

Given the fact that the HF patches are closely associated with the VF patches (Section 
|3]), we quantitatively examine the magnetic fiux of the VF patches associated with each HF 
patch. We estimate the magnetic fiux of the HF patches by considering the cross section of 
the fiux tubes as was done in Section 15.11 We identify the overlapping HF and VF patches, 
and plot the ratio of HF and VF fiuxes in figure [TOl There are cases where an individual HF 
patch overlaps with several VF patches, and vice versa. If the VF magnetic fiux is larger 
than the HF fiux, we take the ratio of the VF fiux to the HF fiux (solid line in figure ITOl) . 
If the HF fiux is larger than the VF fiux, we take the ratio of the HF fiux to the VF fiux 
(dashed line in figure [T0|) . The histograms in figure [10] have the peak at 1, indicating that 
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the magnetic fluxes of associated HF and VF patclies are well balanced. 



6. Intrinsic Magnetic Field Strength 



In this section, we examine the Stokes V amplitude line ratio V6302/V6301 to infer 



intrinsic field strength (jStenflo 



2010l ). Figure [TT] shows the scatter plots of the maximum 
Stokes V amplitudes in the blue lobe of Fe I 630.15 nm and 630.25 nm lines. We measure 
the line ratio of the the pixels inside the orange circle shown in the figure M (c). If the 
intrinsic field strengths are much weaker than 1 kG, the Stokes V is proportional to 
—gdl/dX, where g is the Lande factor ((76302 = 1.667, (^gsoi = 2.50) . In this case, we expect 
the data points to regress to the dashed line in figure [HJ The slope of the dashed line s is 
given by, 

~5'6302<9//<9A|6302 , /^n 

s = ar/aw = ^•^^^ 

"fl'esoiC-' /C'^lesoi 

The value of the slope s = 1.69 is based on the global averages over all the pixels in the 
analyzed subfield of the maximum amplitudes of dl/dX in the blue lobe of Fe I 630.15 nm 
and 630.25 nm lines. 

The data points of the VF pixels inside the HF patches are clearly concentrated along 
the line with the slope of s = 1.69 (figure [11] (a)). Where the VF pixels do not overlap with 
the HF patches, most of data points remain around the dashed line. However, some points 
with larger Stokes V amplitudes start to bend away from this line, and follow the line with 
the slope of 0.7s, shown with the solid line (black contours in figure [11] (b)). This transition 
is due to the differential Zeeman saturation for stronger fields bet ween two Fe I 630.2 lines. 



and the slope of 0.7s indicates the field strength of order of 1 kG (iStenfio 



201( 



Figure [TT] (c) shows the scatter plot for the network fields indicated with yellow 
contours in figure [6] (b). The points are completely away from the line with the slope 
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Fig. 10. — Histograms of magnetic flux ratio when the VF flux is larger than the HF flux 

(solid line) and the HF flux is larger than the VF flux (dashed line). 
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s = 1.69, and regress to the line with the slope of 0.7s (soli d line) a nd 0. 6g for higher 



Stenflol (120101 ) for the strong 



Stokes V signal. The coefficients of 0.6 and 0.7 are given by 
saturated intrinsic field strength. This scatter plot clearly indicates that the network region 
consists of strong kG fields. These strong kG v ertical fields l ocated in th e internetwork 



region would be formed by conyective collapse (jParker 



1978 



Hasan 



1985 



2008; 


Fischer et al. 


2009; 


Stenfio 


2010) 



Nagata et al. 



7. Spatial distribution internetwork magnetic fields on the supergranular scale 

In this section, we qualitatively examine the spatial distributions of the internetwork 
horizontal and vertical magnetic fields on the supergranular scale, and aim to address 
the question of whether the supergranulation affects the distribution of these small-scale 
internetwork fields. 

The area density of the VF patches c?vf(^) at a distance r from the center position of 
the supergranular cell (cross mark in figure |6](d)) is given by 

/ N '^Vf('") / N 

dMr = 9 

where, at a distance between r and r + dr, naii(r) is the total number of pixels, and nyp{r) 
the number of pixels occupied by the VF patches. Here we take dr = 1 pixel (~ 0'.'16). 
When we obtain naii(r) and nYF{r), we do not take into account the network pixels marked 
with yellow contours. Likewise, the area density of HFs di{p{r) is given by 

dHF{r) = (10) 

where nupif') is the total number of the HF pixels marked with red contours at a distance 
between r and r + dr from the center position of the supergranular cell. 

Figure [8] (b) shows the area density dyp^r) and dup^r). Five dashed lines indicate 
r = 4", 9", 15", 22", and 34", and correspond to the orange circle (r = 15") in figure [6] (c) 
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Fig. 11. — Scatter plot of the maximum blue lobe Stokes V amplitudes derived for the 
630.15 nm line vs. the corresponding Stokes V amplitudes for the 630.25 nm line. The 
dashed line with slope s = 1.69 represents the regression relation that would be expected for 
weak fields. Solid lines have the slope of 0.7s. (a) the VF pixels that coincide in position 
with the HF pixels inside the orange circle in figure [6] (c). (b) the VF pixels not overlapping 
with HF pixels inside the orange circle in figure E] (c). (c) network fields shown with yellow 
contours in figure |6] (d). The contours indicate the density of data points. The yellow (red, 
blue and black) contours indicate the area with the density higher than 0.3% (0.1%, 0.05% 
and 0.01%) of the total number of data points. 
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and four circles in figure [6] (d). We find a significant difference in the spatial distributions 
between the HF and VF patches. The area density of the VF patches, dyplr), increases 
with distance. Meanwhile, the area density of the HF patches dn-pif^) does not change from 
the center to the network boundary, and this clearly shows that the HF distribution does 
not have any correlation with the supergranular scale. 



8. Discussion 

8.1. Spatial coincidence of VF and HF patches 

One of the discoveries presented in this paper is the clear positional association of 
the VF and HF patches in the internetwork region. The HF patches are located around 
the boundary of the VF patches as shown in Figure [H] (a), and are connected to the VF 
patches. 96% of the total number of the HF patches has the associated VF patches, 
and 70% of the total area occupied by the HF patches overlaps with the VF patches. 
Alternatively, 49% of the total number of the VF patches has the HF patches. The size of 
the HF patches are apparently smaller compared to that of the VF patches as shown in 
figure [6] (a). Only the small core of the HF patches may be seen due to the lower sensitivity 
of the linear polarization to the magnetic field strength. We quantitatively compare the 
individual VF flux with the flux of the corresponding HF patch, and find that the VF flux 
is close to the HF flux (figure ITOl) . The HF patches so far have been investigated by many 
authors as a new magnetic entity different from the vertical internetwork fields. However, 
they are closely connected, possibly causally related to the well-known VF patches in the 
internetwork regions. The implication of this observational fact will be discussed in the 
subsequent sections. 

As shown in figure [S] (b), the area densities of both the HF and VF patches are 
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essentially constant within r < 15". This is consistent with the observation that the VF and 
HF patches are closely connected. The area density of the HF patches stays the same with 
radius from the center of the supergranule, while that of the VF patches gradually increases 
closer to the supergranular boundary. This behavior at r > 15" is apparently different 
from the correspondence between the HF and VF patches in the pure internetwork region. 
Apparent large-size network fields indicated with yellow contours in figure (b) are not 
included in figure [8] (b) . The histogram of the intrinsic field strengths and the occurrence 
rates of horizontal fields are the same between the quiet Sun and the weak plage region 



when the strong persistent vertical fields (network fields) are removed (jlshikawa fc Tsuneta 



2009aJ). Furthermore, the histograms of the intrinsic field strengths are the same between 



the quiet Sun near the equatorial l imb and the p olar region after removing the vertical 



fields including the network fields (llto et al. 



2010|). The distri butions of the HF and VF 



pat ches in figure M ( b) are consistent with the earlier results by 



and 



Ishikawa fc Tsunetal ( l2009a[ ) 



ItoetaL 



torn . 



8.2. Comparison of longitudinal and transverse apparent flux density 



We estimate and compare the total magnetic flux for both VF and HF patches in the 
internetwork region. The mean flux density for unsigned longitudinal (vertical) component 
is estimated to be 8.3 Mx cm~^ and the mean apparent flux density for transverse 
(horizontal) compon ent is 7.1 M x cm~ ^. We do not find the significant imbalance between 
them as reported by iLites et al.l fcOOSl ). As seen in figure [21 we expect that the photon 
noise contaminates the real signal within Satot or 4(To, making it difficult to distinguish the 
real signals from the noise. Thu s, we exclud e thes e pixels and calculate the mean apparent 

( 120081 ) used all the pixels for the determination 



flux density. On the other hand. 



Lites et al. 



of the mean apparent flux densities in the FOV; their la noise levels are 2.4 Mx cm ^ and 
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41 Mx cm ^ for the longitudinal and transverse apparent flux densities. 



The total unsigned vertical magnetic flux in the internetwork region defined by the 
orange circle in figure |6](c) is estimated to be 3.0x10^^ Mx. In order to accurately estimate 
the total magnetic flux density for the HF patches in the intern etwork region, we con sider 



Ishikawa et al. 



fcoioh . and 



that all the HF patches have the flux tube structure as shown in 
take into account the cross section of the each horizontal flux tube. Indeed, the HF patches 
are considered to be a manifestation of the small-scale bipolar flux tubes as discussed at the 
end of section HI 

The estimated total magnetic flux of the HF patches is 4.7x10^^ Mx, which is 
smaller than that of the VF patches. One reason for this difference may be different 
sensitivity between Stokes V in detecting the longitudinal (VF) component and Stokes 
Q/U in detecting the transverse (HF) component. If emerging loops leave the photosphere, 
there should be more leftover vertical flux, which corresponds to the footpoints of these 
loops, than the horizontal flux. This is both possible and plausible to have more vertical 
flux in addition to the sensitivity issue. Indeed, it is observed that some small-scale 



l oops emerge and reach above the photosphere, leaving 



( iMartmez Gonzalez fc Bellot Rubio 



2009 



Ishikawa et al. 



jehin d the bipolar footpoints 



20101 ). 



8.3. Origin of internetwork fields 



Isobe et al. 



fcoOSh and 



Steiner et al. 



( 2008 



) carried out the numerical simulation of 



the surface layer of the Sun having an initial condition of vertical unipolar magnetic fields. 
They then compared the simulation with observations. Their initialization is similar to 
plage and network regions where vertical fields are dominant. They found that the initially 
uniform magnetic fields are immediately bent by the convective flows, and they eventually 
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become disor dered due to t he co nvective motion and appear as horizontal fields in the 
photosphere. 



Steiner et al 



( 120081 ) also carried out the simulation initialized with uniform 
horizontal magnetic fields, and found that the amount of granular-scale horizontal fields 
exceeds that of the initialization with uniform unipolar magnetic fields. 

We find that the area density of the HF patches does not increase toward the 
supergranular boundary, even though the area density of the VF patches increases toward 
the supergranular boundary. This independence of the HF patches from the amount of 
the VF patches indicates that the seed fields for the HF patches do not come from the 
pre-existing vertical magnetic fields, and our observations appear not to be consiste nt 
with the initializati on with a unipolar vertical magnetic field that 



Isobe et al. 



Steiner et al. 



(12008f l and 



( I2OO8I ) used. Although we cannot rule out other possibilities for the origin 
of the horizontal magnetic fields, it would be conceivable that the granular convective 
motions amplify the field strength of the seed fields up to the equipartition field strength 
corresponding to the granular motion by stretching and folding, and drive the emergence of 
the horizontal magnetic fields from the convection zone. 



8.4. Origin of kG vertical magnetic fields 



2OIOI), we 



Measuring the Stokes V line ratio of the Fe I 630.15/630.25 system (IStenfid 
obtain the intrinsic field strength in an internetwork region. The intrinsic field strength 
inside the HF patches is dominantly weak (< 1 kG), and is in the regime of the weak 
field approximation. Most of the pixels outside the HF patches (but located inside the VF 
patches) also have weak intrinsic field strength, but a significant population of strong (kG) 
fields may also be seen for those regions (figure [TT] (b)). Such signature of kG intrinsic fields 
cannot be found inside the HF patches, as clearly shown in figure [TT] (a). The intrinsic field 
strength in the designated network region is completely above ~1 kG (figure [TT] (c)). 
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We point out that the observed weak intrinsic field strength inside the HF patches 
is in the sub-equipartition regime of the magnetic fields strength. The equipartition field 
strength corresponding to the granular convection is estimated to be 700 G, where the 
density is ~ 10"^ g cm~^ at a depth of 500 km below the surface and the velocity is 
2 km s^^. On the other hand, some pixels outside the HF patches have field strength above 
the equipartition value. 

As shown in figure [9] (d) the HF patches are statistically located inside the granular 
cells, while the VF patches in the internetwork regions prefer to be located in the darker 
region, i.e., closer to the intergranular lanes, though there is a large spread in intensity. 
We conjecture that the magnetic fields below or equal to the equipartition field strength 
are provided by the emergence of granular-scale horizontal fields. When the emerged fiux 
tubes successfully leave the photosphere with left-over footpoints (vertical fiux tube), the 
footpoints are intensified beyond the equipartition field strength (super-equipartition field 
strength of 1-2 kG) due to convective instability. We interpret that the scatter plot of 
figure [m (b) shows such evolutionary process from sub-equipartition to super-equipartition 
fields: outer higher f2— loops have the photospheric footpoints (VF patches) without the 
HF patches, since the top of such loops has already left the photosphere. We speculate that 
such vertical fields corresponding to the VF patches may remain in the sub-equipartition 
regime, or may be subsequently subject to convective collapse, resulting in kG-vertical 
fields. Alternatively, the footpoints of the HF patches may almost immediately become 
strong (super-equipartition) vertical fields due to convective collapse when the HF patches 
are advected by the granular fiows to the intergranular lanes. 



Nagata et al. 



( I2OO8I ) found the observational evidence of the formation of kG fiux 



tube induced by convective collapse using SOT/SP (see also iFischer et al. 



20091); an 



equipartition field strength fiux tube with the preceding transient downfiow of 6 km s ^ 
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goes through the intensification of the field strength to 2 kG. The strong downward flow 
comparable to the local sound speed is observed in the footpoint region of the small-scale 
emerging horizontal fields, suggesting the occurrence of convective collapse at the footpoints 



of horizontal fields ( Ishikawa et al. 



20081 ) ■ However, we point out that there is no direct 
piece of evidence that all the kG vertical magnetic fields in the internetwork region come 
from the emergence of horizontal magnetic fields. Further investigations have to be done. 



8.5. Do horizontal fields concentrate on the supergranular boundary? 



The appearance of the transient horizontal magnetic fields is not spatially random, 
and has an organized cellular structure with a typical scale of 5-10". The cellular 
structures coincide in position with the negati ve divergence of the horizo ntal fiow field. 



i.e., mesogranular boundaries with downfiows (ilshikawa fc Tsuneta 



intern etwork vertical fields, a similar concentration has been observed ( iDommguez Cerdefia 



20101). As for the 



20031 ). Indeed, both the HF and VF patches concentrate in the cellular structure with a 
scale of 5-10" in the snapshot image of figure [6] (a). These structures would correspond to 
the mesogranular cells. A natural question here is whether there is a similar concentration 
of horizontal magnetic fields in the supergranular boundaries. We find that the area density 
of the VF patches increases toward the supergranular boundary, while the area density of 
the HF patches is essentially constant with distance. The HF patches are not concentrated 
on the supergranular boundary. 



The horizontal fiow ve 



300-500 m s-i flShine et al 



o city of the supergranulation at the photosphere is in the range 



20001). Collapsed vertical fields (Section l8.4p can be advected 
by the supergraular fiow fields, and might form the network fields around the supergranular 
boundaries. The abundant vertical fields around the supergranular boundary might be 



originally provided by the emergence of the horizontal fields in the internetwork region. 
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Martini (119901 ) and ISchrijver et al.l (119971 ) pointed out that the main source of the 
network magnetic fields are the ephemeral region and tha t the contribution o f the 
internetwork magnetic fields is not so significant (see also 



Wang &: Zirin 



19871 ). The 



ephemeral region consists of the emerging bipolar magnetic patches with a size less than 
~ 20". We find that the internetwork region is filled with a lot of small-scale emerging 
bipolar magnetic fields, which are not clearly seen with the previous generation telescope. 
The relationship between the ephemeral region and the transient horizontal magnetic fields 
needs to be clarified. 

Here we speculate that vertical magnetic fields in the internetwork fields are advected 
by the supergranular horiz ontal fiow and that t he network fields would be formed. Using 



Hinode/ SOT observation, 



de Wijn et al. 



( l2008l ) find that apparent fiow velocities of 



magnetic elements (patches) in the internetwork region are slightly biased by 0.2 km s~^ 
in a radial direction toward the network boundary. It is clear that the motion of the 
internetwork magnetic patches is affected by the supergranular fiow fields. However, there 
is no direct evidence showing that the network magnetic fields are a consequence of the 
accumulation of the internetwork magnetic fields. In order to confirm that the internetwork 
fields are indeed advected towards the supergranular boundaries forming the network 
fields, further investigation needs to be done with the high spatial and temporal resolution 
observations using Hinode /SOT. 



Ishikawa fc Tsunetal (120101 ) point out that the seed fields of the horizontal fields are 
advected by the mesogranular fiow fields, and granular convective motion amplifies the field 
strength from below the mesogranular equipartition field strength B^q to granular B^q, 
and drives the emergence of the horizontal fields. If the seed fields are advected by the 
mesogranular fiow fields, the seed fields are also likely to be advected by the supergranular 
fiow: we should have seen the concentration of the horizontal fields around the supergranular 
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boundaries. We do not have a clear interpretation for this asymmetry between mesogranule 
and supergranule in terms of the population of the HF patches. This issue would be related 
to a fundamental relationship between the small-scale magnetic fields and the convection 
with different scale sizes. 

9. Summary and Speculation 

In the present paper, we have analyzed the high spatial resolution spectropolarimetric 
data obtained with the SOT/SP to investigate the properties of the internetwork magnetic 
fields. We have paid special attention to the effect of the photon noise in the analyzed data, 
and chosen the valid pixels minimally contaminated by the noise. Therefore, the present 
paper is limited in scope to the field with sufficiently strong polarization signal. An obvious 
penalty of having the stringent noise criteria is that we do not detect weak magnetic fields 
which provide polarization signal comparable to or weaker than the noise level. We do not 
rule out the existence of such weak fields. We refer the pixels with the significant Stokes 
V signal to the VF (vertical field) pixels, and those with the significant Stokes Q or U 
to the HF (horizontal field) pixels. The main observational results with our immediate 
interpretation are summarized as follows: 

1. We find positional association between the vertical and the horizontal magnetic fields 
in the internetwork region. The vertical magnetic fields correspond to the weak 
vertical magnetic fields known for decades to exist. The horizontal magnetic fields 
are investigated in detail with Hinode for the first time. Approximately half of the 
total number of the VF patches accommodates the HF patches partially or entirely, 
and essentially all of the HF patches overlap with the VF patches. These HF patches 
are typically located around the edge of the VF patches. We compare the VF flux 
with the flux of the corresponding HF patch, and flnd that HF flux is close to the 
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VF flux. These observations suggest that the HF patches are physically and causally 
connected to the VF patches in the internetwork region. It has been discussed that the 
transient horizontal magnet i c fields are a m anifestation of the local dynamo process 



( Ishikawa fc Tsuneta 



2009a 



ItoetaL 



20101 ). The origin of the vertical magnetic fields 



in the internetwork region is not yet known, and this positional coincidence may allow 
us to reach a unified understanding of the internetwork fields in general. 

2. Approximately half of the HF patches have bipolar patches (Stokes V signals), and 
the other half have the unipolar Stokes V signal. The unipolar HF patches have 
both a smaller size and smaller horizontal magnetic fiux densities. We suggest that 
only dominant polarity is seen due to the cancellation of the Stokes V signal within 
the resolution element in the unipolar HF patches, and that the HF patches have an 
essentially bipolar structure. Another possibility is that a small portion of the bipolar 
HF patches may be detected as unipolar HF patches. 

3. The mean longitudinal (vertical) magnetic fiux density is estimated to be 8.3 Mx cm~^, 
and the apparent transverse (horizontal) one is 7.1 Mx cm~^. The former iiumber 



is close to the on e derived in the previous studies (e.g.. 



Lites et al 



1[ 

Orozco Suarez et al.il2007b 



20081 ). but the apparent tr ansverse fiu x 



of 5 smaller than the one obtained by 



Lites et al 



density derived here is a factor 



( 120081 ). In our analysis, the mean 



magnetic fiux density of the vertical fields is almost the same as that of the horizontal 
fields. We also estimate the total magnetic fiux of the HF patches, taking into account 
the cross section of the horizontal fiux tubes, and find that the total magnetic fiux of 
the HF patches is a factor of 3 smaller than that of the VF patches. This excess of 
the total magnetic fiux of the VF patches would indicate that the some portion of the 
vertical fields inside the VF patches reach above the photosphere. 



4. Measuring the Stokes V polarization of the Fe I 630.15-630.25 nm line system (line 



ratio method; iStenflol ( 119731 )). we infer the intrinsic field strength in the internetwork 
region. The intrinsic field strengths inside the HF patches (overlapping area of the 
VF and HF patches) are weak (< 1 kG), while the strengths outside the HF patches 
(non-overlapping area) range from weak to strong (kG). This would suggest that 
convective collapse takes place at the footpoint (vertical portion) of the emerging 
horizontal field that is about to leave or that has left the line-forming region in the 
photosphere, and this is considered to be one of the mechanisms to produce strong 
(kG) fields in the internetwork region. The two distinct flux populatio ns represent the 



pre-coUapsed and the post-collapsed vertical flux tubes (IStenflo 



20101) ■ The snapshot 



image analyzed in this paper contains the two evolutionary states. 

5. It is well known that the VF patches concentrate on the mesogranular and 
supergranular boundaries, while we show in this paper that the HF patches 
concentrate only on the mesogranular boundaries of the supergranunlar cell studied 
here. Only the VF patches concentrate around the supergranular boundary. In the 
internetwork region within 15" from the center position of the supergranulation, 
the area densities of both VF and HF patches are constant with radius. The area 
density of only VF patches increases with increasing radius outside 15" radius (closer 
to the network boundary). The area density of HF patches does not increase even 
though there are more VF patches closer to the supergranular boundary, and this 
independence of the amount of the HF patches suggests that the HF patches are not 
due to t he bending of pre- existing vertical magnetic flelds with the granular convective 



motion (jlsobe et al. 



20081 ) 



These observations shed light on the origins of both the vertical and horizontal magnetic 
flelds and their causal relationship, potentially leading to a new picture on the quiet Sun 
magnetism. As a flrst step for a unifled understanding on the quiet Sun magnetism, we here 
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proposed one possible scenario for generation and maintenance of the internetwork vertical 
and horizontal magnetic fields: (1) mesogranular flow fields may advect the seed fields 
with field strength smaller than mesogranular B^q into the mesogranular boundaries. The 
granular convective motion amplifies the field strength of the seed fields up to granular Bcq 
by folding and stretching, and the granular upward motion carries them to the photosphere. 
These are observed as the small-scale horiz o ntal m agnetic fields on the photosphere. This 



process is discussed in 



Ishikawa &: Tsunetal (120101 ). (2) Some horizontal fields reach above 



the photosphere, leaving their footpoints (bipolar magnetic fields). (3) Vertical magnetic 
fields, which are the footpoints of the emerged small-scale loops, are advected toward the 
intergranular lane, and are subject to convective collapse, partially forming strong kG fields. 
(4) Finally, these vertical fields are transported to the supergranular boundaries, and form 
the network magnetic fields. 

The occurrence rate of the e mergence of the horizontal magnetic fields is quite 



high (llshikawa &: Tsuneta 



2009al ). It is likely that such horizontal fields do not 

lere and possibly corona 



return to the photosphere and that they reach the chromosp. 

3 



(iMartmez Gonzalez fc Bellot Rubic 



2009 



Ishikawa et al. 



20101). It is found that 23% of the 



small- scale loops on the photosphere reach abo ve the line forming region of Fe I 630.2 nm 



lines ( iMartmez Gonzalez fc Bellot Rubio 



20091 ). and the total magnetic flux carried by the 



horizontal magnetic fields above the photosphere is e stimated to be t wo orders of magnitu de 



2009 



Ishikawa et al 



20101 ). 



larger than the total magnetic flux of the sunspots (iJin et al.l 
The rest of them would be shredded, thus falling below the threshold of detectability, 
and/or return below the photosphere. It is plausible that the emergence of these ubiquitous 
small-scale loops leaves the vertical magnetic fields on the photosphere, and is the source of 
the quiet Sun magnetic fields with the process suggested above. 
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The idea of this work was inspired by the discussions in the 6th Solar Polarization 
Workshop on 2010 held in Maui. Hinode is a Japanese mission developed and launched by 
ISAS/JAXA, with NAOJ as a domestic partner and NASA and STFC (UK) as international 
partners. It is operated by these agencies in co-operation with ESA and NSC (Norway). 
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